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Abstract 

We study non-equilibrium electronic transport through a quantum dot or an 
impurity weakly coupled to ferromagnetic leads. Based on the rate equation 
formalism we derive the noise spectra for the transport current. We show 
that, due to quantum interference between different spin components of the 
current, the spectrum develops peaks or dips at frequencies corresponding to 
the Zeeman splitting in the quantum dot. A detailed analysis of the spectral 
structure of the current is carried out for noninteracting electrons as well as 
for the regime of Coulomb blockade. 



I. INTRODUCTION 



Resonant transport through a quantum dot (or impurity) has been investigated in nu- 
merous pubhcations. Yet no special attention has been paid to quantum interference effects 
in this process. These effects can generate oscillations in the resonant current through two 
(or more) levels of an impurity. These oscillations are similar to the well-known quantum 
interference effects in the two-slit experiment, and, in turn, produce a peak or a dip in 
the current power spectrum, depending on the relative phase of the two levels carrying the 
current [1]. This feature can be experimentally observed in time-resolved measurements 
of transport currents [2,3]. It has been previously argued that the quantum interference 
effect can explain modulation in the tunneling current at the Larmor frequency in scanning 
tunneling microscope (STM) experiments [4]. 

In this paper we investigate the interference effects in polarized magneto-transport. Con- 
ductance and I-V curves for spin dependent transport through quantum dots has recently 
been studied in several publications [5-9]. Here we study the time dependent properties 
of transport currents. In particular we study the effects related to interference between 
different spin components of the currents. 

These effects can be described schematically as follows. Consider the polarized resonant 
current from the left reservoir (emitter) to the right reservoir (collector) through a single level 
of a quantum dot (impurity) in the presence of an external magnetic field. This field would 
split the resonant level of the dot into a Zeeman doublet. Fig. 1. Let us assume that the 
polarization axis of electrons in the emitter (n) is different from that of the external magnetic 
field (n). Then a spin-polarized electron from the emitter enters into a superposition of 
the "spin- up" and the "spin-down" states of the Zeeman doublet. Fig. 1. As a result, 
the electron wave function in the collector has two components corresponding to different 
energies of the doublet. Yet these components are orthogonal since they correspond to 
different spin components and therefore cannot interfere. If, however, there is an additional 
spin-flip process in the transition between the dot and the collector, the two spin components 



2 



interfere in the collector current. Again, this takes place if the polarization in the collector 
(n') is different from that in the quantum dot (n), Fig. 1. Thus the system operates as a 
two -path interferometer, where the phase difference between the two paths, i.e., through 
the upper and the lower spin states in the dot, contributes to the dynamical (or spectral) 
properties of the collector current. 



Magnetic field 




Fig. 1. Resonant tunneling of a polarized electron through a quantum dot. 
Here the O's denote the tunneling transition amplitudes between the reservoir 
states and the Zeeman doublets (-£^1,2) of the quantum dot. fJ,L.R are the 
chemical potentials in the left and right reservoirs. The unit vectors n, ii and 
n' show the polarization axes in the emitter, quantum dot, and the collector, 
respectively. 

These interference effects can be realized experimentally in a heterostructure with a quan- 
tum dot sandwiched between the two ferromagnetic leads with easy axes different from those 
in the dot. A similar set up can be implemented in other systems, such as self- assembled 
quantum dots [10], ultrasmall particles [11], carbon nanotubes [12] and single molecules [13]. 
These systems are likely to find prominent technological applications, including random ex- 
cess memory and magnetic sensors due to giant magnetoresistance effect [14]. The value 
of Zeeman splitting in the dot is controlled by the external magnetic field. The leads are 
assumed to be thin magnetic films, so that the magnetic field inside the leads is pinned 
parallel to the films and therefore the magnetization in the leads remains unaffected by the 
application of a relatively small external magnetic field. 
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The plan of this paper is as follows. In Sect. II we introduce the model and describe in 
general the rate equation approach for calculations of the resonant current and its noise spec- 
trum. This approach has been obtained directly from the many-body Schrodinger equation 
describing the entire many-particle system [15,16] and allows one to treat the magneto- 
transport through quantum dot systems in the most simple and precise way. In Sect. Ill 
we consider the case non-interacting (of weakly interacting) electrons. We explain in detail 
how the rate equations are used to calculate the polarized current and the noise spectrum 
both for ferromagnetic and non ferromagnetic leads. The results for the time-dependent 
polarized current are compared with the results of the single-particle model, presented in 
the Appendix. We explicitly demonstrate how the polarized current exhibits oscillations due 
to interference effects. In Sect. IV we concentrate on the interacting case. In particular we 
derive the current spectra in the presence of a Coulomb blockade in the dot. We consider 
separately the collector and the emitter currents, as well as the circuit current. In Sect. V 
we summarize our calculations and briefly discuss the potential implications of our results 
on the noise spectroscopy of quantum dots. 



II. MANY-BODY DESCRIPTION. 

Consider the polarized transport of non-interacting electrons through a quantum dot in 
the external magnetic field. Fig. 1. The polarization axis of an electron inside the dot (fi) is 
different from those in the right and left reservoirs (n and n'). The tunneling Hamiltonian 
describing this system can be written as 

l,s d=l,2 r,s' 

+ X ^disCilcid + X ^drs'Cils'ad + H.c \ + Uca\aiala2 , (1) 

\d,l,s d,r,s' I 

where the spin indices, s, s' = ±1/2 are related to different quantization axes (n and n'. 
Fig. 1), and E^gi denote the energy levels in the reservoirs. The Zeeman splitting of the 
dot is denoted by E^^ where d— 1,2. The last two terms describe the tunneling transitions 
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between the reservoirs and the dot states generated by both the tunnehng couphngs Q, and 
the Coulomb repulsion of two electrons inside the dot. 

All parameters of the tunneling Hamiltonian (1) are related to the initial microscopic 
description of the system in the configuration space (x) . For instance, the coupling ^.ais is 
given by the Bardeen formula [17] 

If ^ 
^dis = -TT- / M^) V Xis{x)da , (2) 

where (f)d{x) and Xis{x) are the electron wave functions inside the dot and the reservoir, 
respectively and is a surface inside the potential barrier that separates the dot from the 
left reservoir. Since the spin quantization axes in the dot and in the leads differ from each 
other, the transition matrix elements (O's) in Eq. (2) depend on the relative angles between 
the dot and the lead polarization axes {9l and 9r for left and right leads respectively). The 
simplest form of the couplings that respects SU(2) symmetry is 

nais^^i€£\OL) and Q^rs' ^ ^rd^sl') (Or) , (3) 

where = ±1/2 denotes the electron spin inside the dot. Fig. 1, d^^^'^\9) is spin rotation 
matrix, 

/ cos I sin f \ 

d"'(9)= , (4) 

\ - sm I cos I / 

and fli/r is spin-independent part of the couplings. Neglecting the energy dependence of 
these couplings, Q^^^ = ^l,r, one can relate them to the partial widths (tunneling rates) as 
^L,R = '^t^^\,rPl,r, where pl,r is density of the states in the left (right) reservoir. 

In the case of large bias, \Ei^2 — A*L,i?| ^ ^l,r-i the many-body Coulomb repulsion effects 
in the magneto-transport can be accounted for in the most simple and precise way by using 
the modified Bloch-type equations for the reduced density matrix [15,16]. These equations 
can be derived from the many-body Schrodinger equation by integrating out the reservoir 
states in the hmit of weak or strong Coulomb repulsion, Uc <^ /J'L — Ei or Uc ^ /J'L — Ei, 
without any stochastic or other approximations. In addition these equations are very useful 
for evaluating of the shot-noise power spectrum. 
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In order to apply our method we first redefine the vacuum state, |0), by identifying it 
with the initial state of the entire system. For instance we can identify it with empty dot, 
with the emitter and collector filled up to the chemical potentials Hl,r, respectively. We also 
assume that the electrons in the emitter are polarized along the n-direction. Fig. 1. The 
many-body wave function, describing the entire system can be written in the most general 
way as 



mt)) 



bo{t) + Yl hdis{t)a\ais + hrs>is{t)al^,ais 

d,l,s l,s,r,s' 

^drs'lslsr 



|0), (5) 



where d — {1, 2} denotes a state with one electron in the dot and ls{rs') denote the electron 
level in the emitter (collector). The amplitudes ba{t) of finding the entire system in the 
state "a" are obtained from the Schrodinger equation, idt\'if{t)) = 7i\'if{t)) with the initial 
condition ba{0) — Sa,o 

Let us introduce the (reduced) density matrix aj^jfit), where n, m denote the number 
of electrons arriving the right reservoir with the spin components s' = ±1/2, respectively. 
Figs. 1,2. The lower indices, j, / denote the discrete states of the quantum dot. For instance, 
in the case of non-interacting (or weakly interacting) electrons — {0,1,2,3}, Fig. 2. 
This density matrix, o'™j',"'{t) can be easily constructed from the amplitudes b{t), Eq. (5). 
For example. 



I I 1,1 

CTi2\t) = Ybm/2ma/2it), a','^{t) = Y\bm/2ri/2{t)\\ ... (6) 

I l,r 

The diagonal density matrix elements, crjj™, are the probabilities of finding the system in one 
of the states shown in Fig. 2 and the off-diagonal matrix elements ( "coherencies" ) describe 
a linear superposition of these states. 
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■TkTif ' HTi^ni^ BTbn^ ■r|n^ 

(0) (1) (2) (3) 



Fig. 2. Four available states of the quantum dot. The indices n, m denote 
the number of electrons with the spin components s' = ±1/2 in the right 
reservoir. 

It was demonstrated in Ref. [15,16] that the Schrodinger equation for the entire system, 
idt\^(t)) = H\'^{t)), can be reduced to Bloch-type rate equations describing the reduced 
density- matrix ajj7^{t). This reduction takes place after partial tracing over the reservoir 
states. It becomes exact in the limit of large bias without the explicit use of any Markov-type 
or weak coupling approximations [18]. In the general case these equations are [16] 

\ k k / 

—'Y^V2T^ p{(y jk^k^k'^k' ^j' + Cikji^k^k'^k'-*j) + X/ '^'i'^ P {^k-*j^k'~^j' + ^k^j'^k'^j)o'kk' , 

h h' h h' 

(7) 

(for simplicity we have omitted the indices m and n, which, however, can be easily restored 
from the conservation of the total number of electrons). Here ^Ik^k' denotes the single- 
electron hopping amplitude that generates the k ^ k' transition. We distinguish between 
the amplitudes VL describing single-electron hopping among isolated states and Vl describing 
transitions among isolated and continuum states. The latter can generate transitions be- 
tween the isolated states of the system, but only indirectly, via two consecutive jumps of an 
electron, into and out of the continuum reservoir states (with the density of states p) . These 
transitions are represented by the third and the fourth terms of Eq. (7). The third term 
describes the transitions [k k' j) or {k k' —>■ j'), which cannot change the number 
of electrons (n, m) in the collector. The fourth term describes the transitions {k j and 
k' j') or {k j' and k' j) which increase the number of electrons in the collector by 
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one. These two terms of Eq. (7) are analogues of the "loss" (negative) and the "gain" (posi- 
tive) terms in the classical rate equations, respectively. Yet, the sign of these terms depends 
on the relative sign of the corresponding couplings fl [1]. In our case it is determined by 
the sign of the spin-flip amplitude, Eq. (3). In addition, there is a (permutation) operator, 
7^2 = il, due to anti-commutation of the fermions operators, a| 2 in Eq. (5) (See also [15]). 
The pref actor V2 = whenever the loss or the gain terms in Eq. (7) are generated by a 
two-particle state of the dot. Otherwise 7^2 = 1- 

III. NON-INTERACTING ELECTRONS. 

Consider first the case of no electron repulsion inside the dot, Uc = 0. (In fact, the 
results would be the same ior Uq <^ /J-l — Ei, assuming the couplings Q are independent of 
energy). As in the previous section we choose the initial ("vacuum") state corresponding 
to the polarized electrons in the left reservoir, s = 1/2 (Fig. 1). In this case all four 
configurations shown in Fig. 3 contribute to Eqs. (7). Taking into account that there is no 
direct couphng between the states, £"1,2, i-e. = 0, one obtains the following Bloch-type 
rate equations for the density matrix cr^pit) 

■n,m _ p „n,m , t^W / n-l,m , n,m-l\ , p(2)/ Ji,m-1 , n-l,m\ , p(2) ri,m , n,m 

. Tn(12)/ ?i,m I Ti,m\ -1^(12)/ n— l.m . n— l,m n.m— I n.m— 1\ /o \ 

cTir ^-{rR + 2r?)) a^r + r?VoT + ^sT) - r^'Vrr + ^2T) 

+ rgVaV^''" + rW<73T-' (8b) 

(^22 — - [i- R + ^i- L ) (^22 + L 1^00 +^33 )-^L 1^12 +(^21 ) 

I -nil) ji— l.m . -r\(2) n,m— 1 /o \ 
+ ^r(^33 +^r(^33 (8c) 

■ n.m /oTi I Ti \ I "n(2) n.m , t^{1) n.m , t-i(12) / n.m , n,m\ /oj\ 
(^33 ^ -{^^R + ^L)CT3i +^L(Tii + F^ ^2^ {(^12 +(^21 ) (8d) 

■ n.m (■ I -i-iN n.m -1-1(12)/ n.m , n.m , n.m , n,m\ , t-i(12) / n—l,m n,m—l\ /n \ 

(^12 =-(^e + F)(7l^ -^L {(^06 +0-11 +(^22 +0-33 ) + ^R {(^33 -(^33 ), (8e) 
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where F ^ Fl + Fr and F^l'j, = Fl,rCOs\0l,r/2), F^j, = FL,Rsm''{0L,R/2), F^l% = 
FL,Rsm{6L,R/2) cos(6'l,_r/2) are the partial tunnehng widths of tlic levels £'1,2- 

We now trace the origin of each term in these equations taking as an example Eq. (8b), 
corresponding to j — f — 1 in Eq.(7). The first term in this equation is a "loss" term 
generated by the transitions 1 ^ ^ 1 and 1 ^ 3 ^ 1 in Eq. (7). corresponding to 
the following processes: (a) an electron at the level Ei (Fig. 3(1)) tunnels to the right 
reservoir and back to the same state, with the rate Fr; (b) the same electron tunnels to 
the available continuum states of the left reservoir and back to the level Ei, with the rate 
r^-* — sin^(^^/2). This can proceed only via spin- flip, since the spin-up states in the left 
reservoir are occupied; (c) an electron from occupied states of the left reservoir tunnels to 
the unoccupied level E2 of the dot and then back to the same state of the left reservoir, with 
the rate F^^ 

The second term in Eq. (8b) is a "gain" term generated by the transitions ^ 1, — > 1 
and 3 ^ 1, 3 ^ 1 of an electron from the left reservoir to the level Ei and from the level E2 
to the unoccupied (spin-down) continuum states of the left reservoir. 

The third, "loss", term in Eq. (8b) is generated by the transitions 2 ^ — > 1 and 
2 ^ 3 ^ 1 via the left reservoir. These transitions involve the following processes: (a) 
an electron at the level E2 (Fig. 3(2)) tunnels to an unoccupied, spin-down state of the 
left reservoir, and then makes a spin-flip transition to the state Ei of the dot. The rate of 
this process is (1/2)Fl sin(^^L) cos(^^l/2), as follows from Eq. (7); (b) an electron from one 
of the occupied states of the left reservoir tunnels to the state Ei with the corresponding 
amplitude QlCOs{9l/2). Then an electron with energy E2 tunnels to the vacant state of the 
left reservoir with the spin-flip amplitude — sin(6'i/2). Since this transition proceeds via 
the two-electron state of the dot, the corresponding permutation prefactor, V2 = —1- As 
a result, the rate of this (loss) process is (l/2)Fi sin(^i:). Similar transitions via the right 
reservoir cancel. Indeed the electron from energy level E2 can reach the level Ei by two 
ways: the first through the spin-flip hopping to the right reservoir and then to the level Ei 
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with no spin-flip, and the second, without spin-flip to the right reservoir, and then to the 
level El with the spin flip. These two amphtudes are of the opposite sign. 

The last two terms of Eq. (8b) are "gain" terms generated by the transitions 3 — > 1, 3 — > 1 
of an electron from the state in Fig. 3(3) to the spin-up or spin-down states of the right 
reservoir. The number of electrons in the right reservoir increases by one. 



A. Resonant current in the collector 



Using Eqs. (8) we can easily obtain the spin-up and spin-down currents, Ii/2{t) — 
T,n,miT'Pn,m{t) and I-i/2{t) = lln,m'^Pn,m{t) , whcrc Pn,m{t) = EjS f^ij"" is ^^e probabil- 
ity of finding n electrons with spin up and m electrons with spin down in the right reservoir. 
One finds 



I±l/2{t) 



R 



1 ± COS Or / ,x , 1 T cos Or . . . . sin Or 

7, ^llW 7, (^22{t) + (T33(t) T -^—{ai2{t) + a2l{t)) 



(9) 



where <Jjj'{t) = J2n,m '^Ij^i't)- The latter can be obtained from the following matrix equation 



X{t) + BX{t) = 0, 



(10) 



obtained by the summation of Eqs. (8) over n, m. Here X = {uoo, Cn, C22, ^33, ctu, C21} and 
B is the corresponding 6x6 matrix. 



B = 



( 




-rf - Vr 





r(i2) 


_ra2) X 


r(i) 


Tr + 2rf 





-Vf - Tr 


p(12) 


r(i2) 







Tr + 2T^l^ 


- Tr 


p(12) 


p(12) 





r(2) 


~^ L 




r(i2) 


r(i2) 


p(12) 


p(12) 


p(12) 


p(12) 


ie + r 





[ r(i2) 


p(12) 


p(12) 


r(i2) 

J- L 





-ie + r J 



(11) 



Solving Eqs. (10) and substituting the result into Eqs. (9) we find the following simple 
expressions for the average polarized current: 

TlTr, 



I±l/2{t) 



2V 



;i±cos^LCOseR)(l -e"'*) 

rj^r^rsingLsingfi 

2(e2 + P) 



1 - e"^ ^ cos(et) + e"^ *- sin(et) 



;i2) 
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The same result can be obtained in the framework of a single electron approach, valid for 
the noninteracting case. (See Appendix A.) As expected, the polarized resonant current 
displays damped oscillations. An example of these oscillations in /i/2(^) is shown in Fig. 3. 
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Fig. 3. Spin-up and total resonant currents through the Zeeman doublet as a 
function of time for Ol = Or = 7r/2 and F^^ = = O.le. 

These oscillations, however, disappear in the total collector current, 

m = h/,{t) + /_v2(t) = ^^(1 - e-'') , (13) 
even though electrons in the emitter are polarized. 



B. Noise spectrum of the total collector current 

Now we evaluate the noise-spectrum of the total current, represented by a sum of the 
spin-up and spin-down currents in the final state, Eq. (13). We introduce the density 
matrix crfj,{t) — cr"'^~"'(i), obtained from Eqs. (8), where N denotes the total number 
of electrons which have arrived at the right reservoir by time t. In order to calculate the 
shot-noise spectrum we use the McDonald formula [19] 

S{uj) = 2e^uj / dtsiTi{ujt)—Y,N^PN{t)] , (14) 
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where PN(t) — Yjj=o(^jj(t). One easily finds from Eqs. (8) that 

Y.N'P^{t) = TnY.i^N+l)[a^,{t) + a^,{t) + 2ag{t)] . (15) 

N N 

Substituting Eq. (15) into the McDonald formula, Eq. (14), we finally obtain 

S(u;) = 2e'u;rRlm [Zn(u;) + ^22(0;) + 2^33(0;)] , (16) 
where Z^cu) is a 6- vector, Z = {Zqq, Zu, Z22, ^33, -^12 -^21}, defined as 

%(^') = / E(2^ + (^) exp(^cc;^)d^ . (17) 
One can find Zi^iuS) directly from Eqs. (8) by performing the corresponding summation over 
N . As a result one obtains 

{B - iuI)Z{u) = X + . (18) 

Here B is given by Eq. (11) and / is the unit matrix. The 6- vector X corresponds to the 
stationary solution of Eqs. (10), X ^ X{t ^ 00) and Y{uj) = {Fn + ^22, >33, >33, 0, 0, 0} 
where F (a;) = {loo, ^ii, ^22, 133, ^33, ^12, 121} is given by the equation 

{B - iuI)Y{u) = X (19) 

Using Eq. (18) we calculate the ratio of the shot-noise power spectrum to the Schottky 
noise, S{uj)/2el (Fano factor), where I — I{t ^ 00) = TLTR/Tt, Eq.(13). In particular, the 
result has a simple analytical form for a symmetric dot, = Fj^ = F. We find 

S{u) ^ 2F^+a;^ F V sin^ 

2el 4F2 + a;2 ^ (4F2 + uj^) {AT^ + cj^) ' ^ ^ 

As expected the shot-noise spectrum does not display any peak or dip at frequencies cor- 
responding to the Zeeman splitting, since the interference effects are canceled in the total 
collector current. Yet the noise spectrum depends on the initial polarization of incoming 
electrons {Ol), whereas the total collector current does not (see Eq. (13)). If electrons are 
initially polarized along the magnetic field inside the dot {n) , the Fano factor is the same as 
in the case of resonant tunnehng through a single level [20]. With increasing 9l, however, 
the current flows through both levels of the Zeeman doublet. This leads to an additional 
contribution to the shot noise, described by the second term of Eq. (20). 
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C. Ferromagnetic reservoirs 



Let us consider ferromagnetic reservoirs polarized along n and n' directions, Fig. 1. In 
this case the rate equations (8) have to be modified since there are no available spin-down 
states in the left and right reservoirs. One easily obtains the following rate equations for 
the density matrix where n denotes the number of electron, arriving at the collector 

before time t: 

<o = -rL^o-o + r^^fi + rgV^v^ - r^nHa^,-' + ^^f^) (21a) 
= - (r?^ + r5,^0 <i + r^'^oo - ^(rf ^ - rg^))^^ + a^,) + rgV3T^ (21b) 
^2-2 = - (rg^ + rgO + - ^(rf ^ - rr^)«2 + ^2"i) + rg^<73V^ (2ic) 
<^3"3 = -r«^3"3 + rf + ri^a^, + rf ) «2 + ^2^) (2id) 

^^2 = - (^e + ^r) a^, - T^'^a^, - i(rf ^ - rg^))« + a^"^) + if V^^' (21e) 

Using these equations wc first evaluate the average current, I{t) = h/2{t) given by Eq. (9) 
with (Jjj'it) — J2n^j'j'{'t)- The latter quantities are obtained from a summation of Eqs. (21) 
over n. As a result Eqs. (21) are reduced to the matrix equation (10), where B is the 
corresponding 6x6 matrix of the coefficients of Eqs. (21). Solving this equation we find 
the average current I{t). For instance, in the case of 9l — Or one finds for the stationary 
current, / = -^(oo) — ^l^r/(^l + Fr). One obtains the same expression for the resonant 
tunneling of unpolarized electrons through a single level. 

The time dependence of the average current, I{t), is displayed in Fig. 4 for symmetric 
and asymmetric dots, F^ = Fr = O.le and F^, = e, F^ = O.le, respectively. Comparing 
with Fig. 2 one finds that the oscillations in the average current are more pronounced in the 
case of ferromagnetic reservoirs. This can be anticipated since the corresponding spin-flip 
transitions via the spin-down states of the reservoirs do not exist. We recall that precisely 
these transitions resulted in the cancelation of the interference effects in the previous case. 
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Fig. 4. The polarized resonant current through the Zeeman doublet with 
ferromagnetic reservoirs and Ol = Or = 7r/2. The solid line corresponds to 
Tl =Tr = O.le and the dashed line to = e and Fr = O.le. 



Now we can evaluate the shot-noise spectrum, S{u!) using the McDonald formula. One 
obtains from Eqs. (14) and (21) 

S{u) = 2e''ulm {r^i'> Z^,{u;) + r^R^Z22{ou) + rRZ^^iou) - r^J'^ [^12 (a;) + Z^iicu)]} , (22) 

where Z{u;) is given by Eq. (18) with the matrix B corresponding to Eqs. (21) and 

Y = {Foo,>^ii,?22,0,Fi2,F2i}. Here Foo = cos^ + sin^ ^^22 - ^(1^12 + I21), 

Fn = sin^^Fas, F22 = cos^ ^^33, and F12 = F21 = ^^33, while Yjj, = Yjj,{uj) are 
given by Eq. (19). 

The corresponding Fano factor is shown in Fig. 5 for the same parameters as in Fig. 4. It 
clearly displays a dip at the Zeeman frequency for a symmetric dot. It reflects the damped 
oscillations in the average current, shown in Fig. 4. The dip, however, almost disappears for 
an asymmetric dot with large F^,. 
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Fig. 5. The Fano factor versus lo for a polarized electron current with 
ferromagnetic reservoirs and 9l-,0r = 7r/2. The solid line corresponds to 
Vl = Vji = Q.\e and the dashed line to F^, = e and Vr = O.le. 



IV. COULOMB BLOCKADE 

We now introduce strong Coulomb repulsion inside the dot, Uc ^ jJiL — E^, so that the 
state (3) in Fig. 3 is not available. As a result the corresponding rate equations have an 
even simpler form than those found for non-interacting electrons. Consider again the case 
of ferromagnetic reservoirs, where the quantum interference effects are most pronounced. 
The corresponding rate equations for the case of Coulomb blockade can be obtained from 
Eqs. (21) for non-interacting electrons, by eliminating configurations with two electrons in 
the dot. In the following we consider separately the electron current in the right and in the 
left reservoirs. 

A. Collector current 

The electrical current in the right reservoir and its power spectrum are obtained from 
the following rate equation 
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<7o"o = -r^^o + r«<r^ + rgV^v^ - ^ (^^-^ + a^r') (23a) 
<i = -rSJVl^i + rl^Vs, + + <^2"i) (23b) 

y.(12) 

<^22 = -i i? f^22 + i L f^OO + (<^12 + <^2lj (23cj 

/ r \ r^^^^ 

^12 ^-[ie+ -f) <2 - rf ^<7o"o + -|-(ari + <72"2) (23d) 

Using these equations one finds for the average (polarized) current in the collector 

Ht) = rWan(t) + rgV22(t) - r'^n^auit) + a2,it)] (24) 

where the (Tjj'{t) — J2n,m ^^j^i^) obtained from Eq. (10) for X — {(Too, Cn, (722, cri2, (J21}, 
and S is the 5x5 matrix obtained from the coefficients of Eqs. (23). Solving such a modified 
Eq. (10) for 9l = Or, one finds for the stationary current, 

In = Inioo) = J'^l'l (25) 

This expression shows an asymmetry with respect to the widths and r^^, in contrast with 
the non-interacting case. The reason is that an electron enters the dot from the left reservoir 
with the rate 2rL. However, it leaves it with the rate Tr, since the state with two levels of 
the dot occupied is forbidden. 

The shot-noise power spectrum for the collector current is given by 

Sr{u) = 2e2a;Im {rg^Zii(a;) + r^R^Z22{ou) - T^R^Z^iu;) + Z^iiu)]} . (26) 

Here Zij{u;) are obtained from Eqs. (18), (19), where Y — {1^0,0,0,0,0} and Yqo — 

C0S2 f + Sin^ ^^22 - ^(^12 + Y21). 

The results of our calculations of S{u!) for symmetric and asymmetric quantum dots in 
the case of Coulomb blockade are shown in Fig. 6. 
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5f(0 )/2el 




Fig. 6. The Fano factor versus ui for a polarized collector current with ferro- 
magnetic reservoirs and Coulomb blockade and 9l,0r = 7r/2. The solid line 
corresponds to Fl = Fji = O.le and the dashed line to F^, = e and F^ = O.le. 



B. Emitter current 

We now consider the electric current and its power spectrum in the left reservoir. These 
quantities are determined from the density-matrix a^^, , where p is the number of electrons 
that left the emitter before time t, (the number of holes in the left reservoir). The corre- 
sponding rate equations are similar to Eqs. (23). One finds 

^'oo = -^lct'oo + rJJVfi + rgVf, - r5,^^V?2 + (27a) 

< = -r5,^V?, + pWago"^ + ^{<7{, + afO (27b) 

p(12) 

<Tf2 = -rg^^2 + rgV^o"' + -^{'y'u + ^fi) (27c) 

/ F \ f(^2) 
= - [ie +f)'^i2- ^i'^o^oo' + -^{'^li + <^f2) (27d) 

The average emitter current in the left reservoir is given by /^(t) = ri(Too(i), which differs 
from Eq. (24) describing the collector current, /^(t). Yet, as expected, their stationary values 
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coincide, /l(oo) = Ir{oo), Eq. (25). 

The shot-noise power spectrum of the emitter current is given 

5l (u;) = 2e^a;rilm Zoo (uj) (28) 

instead of Eq. (26) for Sji{u), where Zij{uj) are obtained from Eqs. (18), (19). Yet, Y{uj) = 
{0, cos^ ^,sin^ — ^'"^^^ , — ^'°^^ }yoo('^)) in contrast with the corresponding expression for 
Sr{oj). Even though the expressions for Sl,r{oj) are quite different, one finds that the shot- 
noise power of the emitter current is the same as that in the collector current, Sl{uj) = Sr{lj). 

C. Circuit current 

In general the circuit current is given by Ic{t) — aliit) + piR{t), where the coefficients 
a, P with a + P — 1 depend on the junction capacities [21]. Using charge conservation, 
II — Ir + Q, where Q is charge in the dot, one finds 

Ic{t)Ic{0) = alL{t)lL{0) + (3lR{t)lR{0) - aPQ{t)Q{0) . (29) 

Using this relation one finds a simple expression for the noise spectrum of the circuit current 
[4,22] 

Scico) = aSLico) + PSr{u) - aPco^SQiu) . (30) 

where Sq{u!) is Fourier transform of the charge correlation function. This quantity can 
be obtained straightforwardly from the matrix equation (19), where X is the 5- vector 
{0, an(oo), ^22(00), 0, 0}. Then Sq{uj) = 4Re[Fn(c^) + 1^22(c^)]. 

The results of our calculations of Sc{u;) lor a — (3 — 1/2 are shown in Fig. 7. 
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Fig. 7. The Fano factor for the circuit (a = /? = 1/2) versus a; for a polar- 
ized electron current with ferromagnetic reservoirs and Coulomb blockade and 
Ol^Or = 7r/2. The solid line corresponds to Tl = Tr = O.le and the dashed 
line to r^, = e and Vr = O.le. 

One finds from Figs. 6 and 7 that the Coulomb blockade modifies the current spectrum 
very drastically with respect to the non- interacting case, Fig. 5. 

V. CONCLUSIONS 

In this paper, we study the interference effects in magneto-transport through Zeeman 
split levels of quantum dots or impurities. We concentrated on the time-dependent prop- 
erties and the power spectrum of the electric current by applying the new approach using 
quantum rate equations, which is mostly suitable for this type of problems. We explicitly 
demonstrated that our method produces the same results as a single electron approach, 
widely used for a description of non- interacting electron transport. Yet the quantum rate 
equations method is valid also for the case of interacting electrons and accounts for the 
Coulomb blockade in the most simple and precise way. 

Our results indicate that the Coulomb blockade plays an important role in the spectral 
properties of the transport current. First of all, in the presence of Coulomb blockade the 
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signal-to-noise ratio significantly is amplified, as one can observe from the results in the 
previous sections. This is probably a consequence of the prohibition of double occupation 
of the resonant level in the quantum dot. Indeed, when two electrons in the dot are present, 
the interference effects are suppressed due to the "randomization" of the relative phase. 
Interestingly, the dip in the noise spectrum for the noninteracting electrons is replaced by a 
peak as result of Coulomb interaction. Clearly the Coulomb interaction modifies the phase 
of the electrons tunnehng trough the dot, which "flips" the spectral feature in the noise. 
The details of this very interesting phenomena must be studied in the future. 

We emphasize that the coherent oscillations in the current can be observed only for po- 
larized current and that oscillations disappear for unpolarized current. This is different from 
the resonant transport through two orbital levels of a quantum dot or impurity, where the 
quantum interference effects can be observed even in unpolarized case. Therefore it is most 
natural to use ferromagnetic leads for observation and utilization of quantum interference 
effect in the magneto-transport. Thus our calculations were mostly concentrated on this 
case. Our results show explicitly the appearance of peak or dip at a frequency near the 
Zeemann splitting frequency (Larmour frequency). We believe that this phenomenon can 
be useful for analyzing the noise spectroscopy of quantum dots or impurities. Indeed, the 
Zeeman splitting of a localized quantum dot orbital must be sensitive to local magnetic 
fields, and therefore one can hope that such coherent effect, if observed experimentally, may 
allow for detection of the local hyperfine structure of the dot/impurity. This, however, must 
be a subject of a separate investigation. 
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APPENDIX A: SINGLE-ELECTRON DESCRIPTION 

In the case of non-interacting electrons one can compare our results with those obtained 
using a single-electron approach. Although the latter is widely used in the literature, it is 
usually restricted to the time-independent (stationary) case. Here we present an extension 
of the single-electron approach for the non-stationary case. This would allow us to evaluate 
the time-dependent resonant current. Fig. 3, and to compare the results with those obtained 
from Eqs. (7). 

Let us consider a system consisting of the reservoirs and the quantum dot filled with 
only a single electron. We assume that this electron is initially in the left reservoir (emitter) 
at the level Ej with the spin polarized along the n-direction. Fig. 1. The electron motion is 
described by a wave function which can be written in the most general way as [23]: 



l,s d=l,2 r,s' 



|0) , (Al) 



where ha{t) is the amplitude of finding the electron in the state a given by a corresponding 
creation operator. These amplitudes are obtained from the Schrodinger equation |^(i)), 
with the initial conditions his{0) — 5//5s,i/2 and bd(0) = brs(O) = 0. It is useful to use the 
Laplace transform, h{E) = /q°° b{t) exp{iEt)dt. In this case the time- dependent Schrodinger 
equation for the amplitudes b{E) becomes the following system of linear algebraic equations 



{E - Eis)hUE) -QiYl €,SiOL)hAE) = iSi,T5s,i/2 (A2a) 

d'=l,2 

{E - eA{E) - Y: niS^J,?{er.)bUE) - ^ 0,di^jJV«)b,,,(£^) = O (A2b) 

{E - Ers')brs'{E) -Qr d^'j (^i?) Brf^^) = • (A2c) 

d'=l,2 

Substituting bis and b^.^' from Eqs. (A2a), (A2c) into Eq. (A2b) and replacing the sums on 
I and r by the integrals, we obtain 
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E-Ei+ I bi (E) = I— (A3a) 

E-E2 + 1 — ^ — j "^'(^^ " ~'~E^r^~;7' ^ ^ 

Note that the amphtudes, bi(£') and h2{E), are decoupled in Eqs. (A3) although the corre- 
sponding states are connected via the continuum. The reason is that the spin-flip couplings 
of the dot with the reservoirs are of the opposite sign for the spin-up and the spin-down 
states of the dot {Ei and E2 in Fig. 1). However, for the general case of resonant tunneling 
through two levels, the corresponding amplitudes are coupled via the interaction through 
continuum [1]. 

Using the inverse Laplace transform bi,2(^) = / bi,2(-E') e:K.p{—iEt)dE / (27r), we obtain for 
the amplitudes, bi,2(t), for finding the electron inside the dot 

^ ^LCOs{ej2) / _ ^A4a) 

b,(t) = _^}L^ihm_ f-iE,t _ ^-iE,t-lt\ (A4b) 

El — E2 + i-2 ^ 

where F = F^ -|- F^j. The probabihty amphtude of finding the electron inside the collector 

is hrs'{t) = J brs'{E) exp{—iEt)dE / {2n) , where hr,s'{E) is given by Eq. (A2c) 

~brAE) -^^j:diZ\(^R)UE) . (A5) 



The above equations determine the motion of a single electron placed initially in the 
emitter. In order to obtain the polarized current, Is'{t), in the single-electron model one has 
to sum over all initially occupied states Ej of the emitter and over all available states Er of 
the collector. Thus Ig' — dNs'{t)/dt, where Nsi{t) — J^lr |brs'(^)P is the average number of 
electrons with spin-up and spin-down (s' = ±1/2), accumulated in the collector by the time 
t. Using the inverse Laplace transform and replacing J2ir ~^ I PLPRdE^dEn we obtain 



AT,, 



it) = J PLPRdELdEn j ^^hrAE)b;AE')e^^'''-''^' (A6) 



Substituting Eq. (A5) into Eq. (A6) and integrating over Ers' one obtains for the polarized 
current 

22 



h/2{t) = Tr / plcIEl] cos(V2)bi(0 - sm{9R/2)h2{t)\^ (A7a) 
/-i/2(i) = r« / pLdEL\ sm{9R/2)hi{t) + cos{9R/2)h2{t)\^ (A7b) 

where the amphtudes bi,2(i) are given by Eqs. (A4). Note that these amphtudes in the 
stationary hmit, bi_2(^ oo), are the transmission amphtudes describing the resonance 
tunnehng through the levels £^1,2 [7]. Thus Eqs. (A7) represent a generalization of the 
Landauer formula for the time-dependent case. 

For large bias, /xl — /i/j ^ F, the integration over El in Eqs. (A7) can be performed 
analytically using Eqs. (A4) for the amplitudes bi^2(^)- As a result we finally arrive at 
Eq. (12) obtained from Eq. (7) for the case of non-interacting electrons. This agreement with 
the case of non-interacting electrons is quite remarkable since our rate equations dealing with 
many-electron states are very different from those obtained in the single electron framework. 
Yet, this is not surprising since in the case of non- interacting electrons the single electron 
description is valid. In fact, Eqs. (A3) can be mapped to Eq. (10) using |bj(t)p = crji(t) -|- 
0-33 (t), where i = 1,2 and hi{t)h*2{t) = auit). 
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